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START: Signal received at \ 
antenna (dE/dt, dB/dt) J 




NO 



Digitize at 20 
MHz at 12 bits 



Mark threshold 
crossing time 
using GPS 



Sum every 4 
samples (5 MHz at 
14 bits) 



IOC 



Store threshold 
crossing time to a 
threshold crossing 
timestamp buffer 
(FIFO) 



( H lo 



Store 5 -MHz 
sample to 1 st 
buffer (FIFO) 
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From 170 
Fig 10(a) 



Transfer data FIFO 
sample to circular 
buffer for DSP 



Get oldest threshold 
crossing timestamp 
from timestamp 
buffer 





Determine onset 
time of pulse by 
looking backward 
until samples reach 
or cross noise Level 



Proceed to Pulse 
Classification 



I 



Add current sample 
to digital numerical 
integration for field 
amplitude 



To 242 
Fig 10(e) 



Initialize digital 
numerical 
integration to 
obtain field 
amplitude 



To 192 
Fig 10(c) 
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To 1 Tl 
Fie 10(b) 



From 190 
Fig 10(b) 



1 



Record integrated (field) 
signal amplitude for this 
sample 



Provide smoothed field 
derivative signal using 
surrounding samples 



Find any peak or trough 

in 




Determine peak time 
and amplitude using 
smoothed derivative 
and field signals 



Increment counter 
for number of 
peaks in pulse 



2oi 

Find all zero crossings 




Record zero- 
crossing time as 
time of previous 
sample 



r*. We) 



To 206 
Fie 10(d) 



lightning Detection and Data Acquisition System" 
Martin Muiphy et aL; Serial No. 10/075,055 



From 204 
Fig 10(c) 



Increment counter 
for number of zero 
crossings 



Continue by reading next 
data sample into DSP CPU 




To 210 
Fig 10(b) 
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From 184 
Fie 10(b) 
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Pulse Classification 




Assume pulse is 
due to cloud 
discharge 



Is pulse likely to be due to noise? 



Discard pulse 




35^ 



Clear counters for 
peaks and zero 
crossings 



NO 



To 174 
Fig 10(b) 



Opposite-polarity 
peak time test 
result false 




Test opposite-polarity peak time difference 



Opposite-polarity 
peak time test 
result true 
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To 260 
Fig 10(f) 
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From 256 
Fig 10(e) 



Test bi-polar amplitude ratio 




NO 



Classify pulse as 
due to a leader 
(sub-class of IC) 



j — 16 8 



Classify pulse as 
bipolar IC (sub- 
class of IC) 



Classify pulse as 
due to a CG return 
stroke 



Proceed to save pulse information 



I 



To 284 
Fig 10(h) 



NO (option I) 



Classify pulse as 
due to distant CG 
return stroke 

V— 270 




Discard pulse 



Test first peak-to-zero time ^ 
YES 



/ '^2 



Reclassify pulse as 
due to a IC 
discharge 




To 174 
Fig 10(b) 
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From 258,272,278 
Figs 10(f), 10(g) 




Save Pulse Information 




DEFINE PULSE: 
1 st & 2nd peak times 
Onset time 

Threshold crossing time 
First risetime 
1st & 2nd peak amplitudes 
1 st peak-to-zero time 
Opposite polarity peak 
time and amplitude 
Field derivative at first 
zero crossing 
Pulse classification 

$ 



Store pulse 
information in 
Results Buffer 





Continue by reading next 
data sample into DSP CPU 



Proceed to Pulse Train Processing 



To 412 
Fig 14 
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To 174 
Fig 10(b) 
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From 262,271 
Fig 10(f) 




due to IC 
discharge 



Proceed to save pulse information 
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To 284 
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From 240, 290 
Fig 10(h) 




Transfer 2 n 
newest pulse to 
host PC 
individually 
with fiill 
feature set 



Remove 2* 1 
newest pulse 
from results 
buffer 



•<Uo 



Transfer pulse 
train to host 
PC interface 



Optional 
PTK test 
(FIG. 15) 



Initiate new 
pulse train 
using 2 
newest pulse 



Remo\e 2 nd 
newest pulse from 
results buffer 



Zero pulse 
train info/Clear 
train buffer 




Flag newest 
pulse as largest 
in train 



Move newest pulse 
in results butter to 
position of 2 nd 
newest; await next 
pulse (up to max 
wait time) 



To 412 
Fig 14 
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Return to 
processing 
where PTK test 
occurred (d. 5) 



To 436 
Fig 14 



Transmit pulse 
individually with 
full feature set 



Discard pulse 



train 



Remove 
newest pulse 
from results 
buffer 



Remove 2 
newest pulse 
from results 
buffer 



Zero pulse 
train info/Clear 
pulse train 
buffer 



To 412 
Fig 14 



"Pre-Trigger Kibosh " test 
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START Pulse Train Compression 



Clear 
amplitude 
fraction 
indicator bit 



Find position of 
largest amplitude 
pulse in train (already 
time-sorted) 



-2- 5M 



Save timestatnp to the 
second (32 bits) and 50- 
nsec counts (25 bits) of 
largest pulse in data 
record. 



SlO 




Set amplitude 
fraction 
indicator bit 



Save full amplitude 
(14 bits) of largest 
pulse in data record. 
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From 518, 520 
Fig 20(a) 



Encode 

classification of 
largest pulse (1 = 
Cloud. 0 = CG) 



Set loop counter to 
position of largest 
pulse minus one. 
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From 540 
Fig 20(c) 




5^ 



Time difference = Time of 
pulse indicated by loop 
counter minus time of 
pulse after it. 



Encode time difference 
(10 of microseconds in 
signed 13-bit value) 



5^S 



53o 




534 



Amp. fraction = 1000 * 
ratio of amplitude of 
current pulse to amplitude 
of largest Dulse 



5%, 




To 538 
Fig 20(c) 



To 542 
Fig 20(c) 
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From 536 
Fig 20(b) 



From 526 
Fig 20(b) 



Encode 

classification of 
current pulse (1 = 
Cloud, 0 = CG) 



To 526 
Fig 20(b) 



if!! 



53S 



Decrement loop 
counter 



Set loop counter to 
position of largest 
pulse plus one 




Time difference = Time of 
pulse indicated by loop 
counter - Time of pulse 
before it. 



Encode time difference 
(lOths of microseconds in 
signed 13 -bit value) 



To 550 
Fig 20(d) 



To 560 
Fig 20(d) 
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From 598 
Fig 20(g) 



From 544 
Fig 20(c) 









Use amplitude 






v. fractions? 








r 




Amp. fraction = 1000 




* ratio of amp. of 




current pulse to amp. 




of largest pulse 




r 




Encode amplitude 




fraction (10 bits) 
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ft! 
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Encode pulse 






classification (1 = 






Cloud, 0 = CG) 








To 544 




Increment loop 


Fig 20(c) 


< 


counter 



DONE. Transmit 
encoded data to 
central analyzer. 



r 
v. 



END Compression 
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If! 



START Decompress Pulse 
Train 



Use communication protocol 
header info to determine total 
size of train data record 



Determine time of largest 
pulse 



57* 



574 



578 



5*6 




>YES 



# pulses = (record length - 
1 1) + 2 + 1 {rec. length in 
bytes} 



# pulses = (record length - 
11) -5- 3 + 1 {rec, length in 
bytes} 



Decode amplitude of largest 
pulse 



Decode classification of 


largest pulse (1 = Cloud; 0 - 


CG) 














To: 588 






Fig. 21(b) 
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From 586 
Fig 21(a) 



Set loop counter to 
zero 



To 614 
Fig 21(c) 



YES 



Decode time 
difference 



Add time diffto 
previous pulse time to 
get this pulse's time 



Add timediff. to 
largest pulse time to 
get this pulse's time 



Decode pulse 
classification (1 = 
Cloud; 0 = CG) 



Fig 21(b) 



To 606 
Fig 21(c) 



5^0 




From 612 
Fig 21(c) 




Add time diff. to 
largest pulse time to 
get this pulse's time 



x: 
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From 604 
Fig 21(b) 




largest pulse amp. to 
determine pulse amp. 



Increment loop 
counter 



To 590 
Fig 21(b) 



414 



END: Decompress pulse train 




Fig 21(c) 
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